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Abstract 

Acousto-elasticity is concerned with the propagation of small-am- 
plitude waves in deformed solids. Results previously established for 
the incremental elastodynamics of exact non-linear elasticity are use- 
ful for the determination of third- and fourth-order elastic constants, 
especially in the case of incompressible isotropic soft solids, where the 
expressions are particularly simple. Specifically, it is simply a matter 
of expanding the expression for pv 2 , where p is the mass density and 
v the wave speed, in terms of the elongation e of a block subject to a 
uniaxial tension. The analysis shows that in the resulting expression: 
pv 2 = a + be + ce 2 , say, a depends linearly on p; b on p and A; and c 
on p, A, and D, the respective second-, third, and fourth-order con- 
stants of incompressible elasticity, for bulk shear waves and for surface 
waves. 
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1 INTRODUCTION 



Recent years have witnessed a surge of interest in the acoustics of incompress- 
ible soft solids, with the long-term goal of mastering all aspects of ultrasonic 
wave propagation in biological soft tissues. In the wake of a 2004 article by 
Hamilton et al. [T], establishing that the general strain energy W of fourth- 
order incompressible elasticity involves only three elastic constants, 

W = n\x (E 2 ) + ^ tr (E 3 ) + D (tr(£ 2 )) 2 , (1) 

where E is the Green strain tensor and /i, A, and D are second, third-, and 
fourth-order elasticity constants, respectively, at least 16 articles (2j El El El 
El [3 El [91 HUl HH H2l [13l HH HSl HZ] have studied the dynamics of those 
solids. 

Although it is of course important to make progress in this area, because 
of obvious repercussions for elastography techniques, it is also important to 
acknowledge, and to use, related and previously established results, such as 
those coming from the literature on incremental (small-on-large) waves in ho- 
mogeneously deformed materials, because they give direct access to a simple 
way to evaluate /x, A, and D. In fact, the acousto-elasticity of non-linearly 
elastic incompressible solids has a long and strong history, dating back to 
the seminal works of R.S. Rivlin, M.A. Hayes, M.A. Biot, and many others. 
For instance, Rivlin and Saunders|18] find in 1951 that only two constants 
are needed to encompass third-order incompressible elasticity (where W is 
given by Eq.(JTJ written at D — 0), and Ogden[T9] shows, as early as 1974, 
that only three constants are required at fourth-order, as in, and preceding, 
Eq.©. 

The present paper shows that acoustic plane waves can be used to measure 
the third- and fourth-order elastic constants of incompressible soft solids. 
This goal can be reached using non-linear waves, but these are difficult to 
generate and to detect experimentally in solids. Moreover, although it might 
be possible to observe them in liquid-like solids p2] (such as gels, phantoms, 
agar, etc.), it seems harder to imagine that they could travel in biological 
soft tissues without causing damage at the cellular level, by virtue of them 
being of large amplitude (see Mironov et al. [20] for experimental evidence of 
such damage). On the other hand, linear, or rather linearized, acoustic waves 
can also give the non-linear elastic constants, through their coupling with a 
pre-deformation: this is the so-called acousto- elastic effect. 
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An abundance of results exists in the literature on incremental dynamics 
in exact non-linear elasticity, usually giving the speed of a small-amplitude 
wave propagating in any incompressible solid, subject to any homogeneous 
pre-deformation. It is thus a simple matter to specialize these results to the 
strain energy density of Eq.(JTJ, and to a uni-axial tension yielding a stretch 
A = 1 + e, say, where e is the measurable elongation. 

When an isotropic soft solid is subjected to a uni-axial tension (e > 0) 
or compression (e < 0), it exhibits strain-induced anisotropy. The axes of 
anisotropy are the so-called principal axes, aligned with the principal axes of 
the Cauchy-Green strain ellipsoid. The squared wave speed of an infinitesimal 
wave traveling in an incompressible solid subject to a small-magnitude uni- 
axial pre-deformation depends (at first) linearly on the elongation e, and 
the coefficients of linear interpolation are linear combinations of p and A. 
Conversely, measuring the wave speed gives access to the values of p and A. 
In order to access D, the fourth-order constant of ([1]), the characteristics of 
a small-amplitude wave must be coupled to a small-but-finite uni-axial pre- 
deformation, i.e. yielding v 2 as a quadratic function of e. In other words, this 
paper aims at finding explicitly the coefficients a, b, and c in the expansion 

p V 2 = a + be + ce 2 , (2) 

where p is the mass density. For bulk homogeneous waves (Section III), 
we find of course that a = p, the initial shear modulus, giving the linear 
wave speed \fpfp~. Similarly, we find that a = 0.9126/x for surface waves 
(Section IV), consistent with Lord Rayleigh's result [21] in linear isotropic 
incompressible solids. For both types of waves, we find that b depends linearly 
on p and A, and that c depends linearly on p, A, and D. 

2 PRE-STRESS, PRE-DEFORMATION, 
ELASTIC MODULI 

Take a parallepipedic sample of a soft, isotropic, incompressible solid and sub- 
ject it to a uni-axial tension or compression, and call A the resulting stretch in 
the direction of elongation. It is a simple and common exercise of non-linear 
(exact) elasticity to show that the corresponding principal stretches are: 

_i _i 
Ai = A, A 2 = A 2, A 3 = A 2, (3) 
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showing that the deformation is equi-biaxial. We call X\ the axis along the 
direction of elongation, and X2, £3, two orthogonal axes in the plane normal 
to that direction, such that the faces of the deformed block are parallel to 
the (xi,X2)-, (x2,xs)-, and (X3, xi)-planes. 

With respect to incremental wave propagation, some quantities, 7^, /3y 
(i, j = 1, 2, 3), play a privileged role, akin to that of elastic moduli. They are 
defined in general by [221 ES] 

7« = (W« - XjWjW/M - A 2 ) ? 7ii, 

Ai = I (Af^n + AJW^) - AjAyWjj - (AjWj - A i W J )A i A i /(A? - A|) = 

(4) 

when A, 7^ Aj, and by 

Tii = I (A?W« - XiXjWij + AiWi) 7^ 7ii , 

N = \ (Ajir /; - AjAjWjv,- - Wi) = (5) 

when Aj = Aj. Here Wi = dW/dXi, Wij = d 2 W/d\idXj, and z 7^ j (also, 
there are no sums on repeated indices). We recall that the eigenvalues Ei of 
the Green strain tensor are related to the Aj by Ei — (Af — l)/2. It is thus 
easy to write W of Eq. ([TJ in terms of the A.; and then to relate the results to 
Eq.((3]). Finally, the results can easily be expanded in terms of the elongation 
e = A — 1. Here we find that 

712 = 7i3 = + (3/i + A/ A) e + (5/x + 7A/4 + 3D) e 2 + C(e 3 ), 

721 = 73i = /i + (A/A)e + (2/i + A + 3D)e 2 + G(e 3 ), 

012 = As = /i + (3/i/2 + A/4) e + (25/i/2 + 47A/8 + 12D) e 2 + C(e 3 ), 

723 = 732 = £32 = fi ~ (3// + A/2) e + (5/i + 7A/4 + 3D) e 2 + C(e 3 ), (6) 

and we do not need to specify the higher-order terms for our purpose, which 
is to establish relationships in the form of Eq.(l2l. Notice that 

2/3i2 - 712 - 721 = 2/3i3 - 713 - 73i = 12/i^e 2 + C(e 3 ), 

3 / A/2 + D\ . . 
where v = - M + — - J (7) 

is the coefficient of non-linearity for linearly polarized plane shear waves, 
obtained by Hamilton et al.pQ 



4 




Figure 1: Deformation of a unit cuboid in an initially isotropic, incompress- 
ible solid, elongated by the amount e = A — 1 (in this sketch, the cube is 
in tension and e > 0). The shear bulk waves under study propagate in ei- 
ther of the three principal planes (here, the (xi, X2)-plane), in the direction 
of the unit vector n, at an angle 9 with a principal direction (here, the x±- 
direction). One wave is polarized in the principal plane along the unit vector 
a orthogonal to n, the other wave is polarized along a x n, the normal to 
the principal plane. 

3 BODY WAVES 

In a deformed incompressible solid, two infinitesimal transverse (shear) ho- 
mogeneous plane waves may propagate (but no longitudinal homogeneous 
plane wave), which are linearly polarized along mutually orthogonal direc- 
tions. Here we consider the propagation of homogeneous plane waves in the 
(xi,X2)-, {x2,xz)-, and (xs, xi)-planes, the so-called principal planes. 

First consider that the waves travel in the (xi, X2)-plane. In that plane, 
call n the unit vector in the direction of propagation and call a the unit vector 
orthogonal to n, see Figure [TJ Based on previous studies (see, for example, 
Ogden[23]), we find that one of the two shear bulk waves is polarized along 
a, and travels with speed v a , given in terms of the elastic moduli by the 
following secular equation, 

pv 2 a = (7i2 + 72i - 2/3 12 ) cos 4 9 + 2{(3 12 - 721 ) cos 2 e + (8) 
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in general, where 9 is the angle between n, the direction of propagation 
and x\, the direction of the uni-axial pre-deformation (both the angle and 
the phase speed are evaluated in the deformed configuration). Using the 
expressions Eq.® for 712, 0i2, and 721, we arrive at a relationship in the 
form of Eq. ([2]) , where 

a = fi, b = 3/i cos 2 9 + A/4, 

c = 2fi + A + 3D + (21/i + 39A/4 + 18D) cos 2 9 - Yl\w cos 4 9. (9) 

The relationship provides an obvious way to determine /x, A, and D exper- 
imentally, either by fixing e and varying 9, or vice-versa, see the example 
presented in Figj2j The other bulk shear wave is polarized along £3 (i.e. 
along a x n), and its wave speed v\, is given by 

pVb = 7i3 cos 2 9 + 723 sin 2 9, (10) 

in general or, using the expressions Eq.([6]) for 713 and 723, also by a relation- 
ship in the form of Eq.©, where now 

a = /i, b = 3(/x + A/8) cos(2#) - A/8, c = 5/x + 7A/4 + 3D. (11) 

Notice here how the quadratic coefficient in the pv 2 -e relation is independent 
of the angle of propagation. Note also that in the case of principal wave 
propagation along x\, the direction of elongation, 9 = and the two shear 
wave speeds coincide: 

pv 2 a = pv 2 b = n + (3/i + A/4) e + (5/i + 7 A/4 + 3D) e 2 , (12) 

showing that x± is an acoustic axis, along which circularly-polarized waves 
may propagate. 

Next, consider that the waves travel in the (x\, X3)-plane. Again, call 
9 the angle between the direction of propagation and x±, the direction of 
the uni-axial pre-deformation. Then we find that, owing to the equalities in 
Eq.flS]), the secular equations are exactly the same as for propagation in the 
(x\, X2)-plane. This is consistent with the transversely isotropic character of 
the strain-induced anisotropy. 

Finally, consider waves traveling in the (X2, X3)-plane. We call 9 the angle 
between the direction of propagation and X2, so that n = [0, cos 9, — sin9] T 
and a = [0, sin 9, cos 9] T . Then the secular equation for the shear wave 
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polarized along a is Eq.([5|) where the indices 12 and 21 are replaced by 23 
and 32, respectively. However, because 723 = 732 = /?32, the dependence on 9 
vanishes, and the wave propagates in an isotropic manner in that plane. Its 
speed is given by 

pv 2 a = 723 = fi - (3/i + A/2) e + (5/i + 7A/4 + 3D) e 2 . (13) 

Similarly, the secular equation for the shear wave polarized along n x a, i.e. 
along x\, the direction of elongation, is Eq. ffTO]) . where the indices 13 and 
23 are replaced by 21 and 31, respectively. However, because 721 = 731, the 
dependence on 9 vanishes and that wave also propagates in an isotropic man- 
ner in the (X2, £3)-plane. This observation is consistent with the transversely 
isotropic character of the induced anisotropy. In that latter case, the speed 
is given by 

pv 2 b = 721 = /x + (-A/4)e + (2/i + A + 3D) e 2 . (14) 

Now we make the link with the theoretical results of Gennisson et al.[S], 
who derived the linear acousto-elastic dependence of the squared wave speed 
on the uni-axial stress a (say), in the cases of principal wave propagation. 
To establish this connection, it suffices to keep the linear part of the pv 2 -e 
relations Eqs. (fT2l . (1151) . and (TT4"]) . and to recall that for incompressible solids 
e = — cr/(3/x) + . . ., to recover the following expansions [S], 

^ = "-'( 1 + i)' + (15) 

respectively. 

We conclude this section with an example taken from experimental in- 
vestigations. Renier et al. (Sj use a small elongation (compression) to deduce 
a linear pv 2 -a relationship from experimental data on Agar-Gelatin based 
phantoms: they conclude that for their 5% Gelatin sample, fi ~ 6.6 kPa and 
A ~ —37.7 kPa. Then they use finite-amplitude shear waves to estimate the 
coefficient of non-linearity: they find v ~ 3.5 (see also the combination of 
Ref.[3] and Ref.[T2].) Eqs.(18])- (jl4p provide an alternative mean of determining 
v for that sample, using infinitesimal^ non-destructive, shear waves instead 
of finite-amplitude shear waves. As seen clearly on FigJH the (theoretical) 
linear part of the pv 2 -e is confined to ±1% stretches, approximatively, and 
it suffices to elongate the block by a few percent to perceive the influence of 
the fourth-order elasticity constant. 
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Figure 2: Shear bulk waves propagating in a uni-axially deformed incom- 
pressible block: pv 2 (in kPa) against the elongation e in the case of waves 
traveling (i) in the Xi-direction (then, all shear waves travel with speed given 
by Eq.(12), irrespective of the direction of polarization); (ii) in any direction 
in the (x\, X2)-plane or in the (xi,x-s) -plane, with transverse polarization in 
that plane, and speed given by Eq.(13); (in) in any direction in the (xi,x 2 )- 
plane or in the (x\, X3)-plane, with polarization normal to that plane, and 
speed given by Eq.(14). The thin straight lines correspond to the linear part 
of the pv 2 -e curves. Here the initial shear modulus is \x = 6.6 kPa, the third- 
order Landau constant is A = —37.7 kPa, and the coefficient of non-linearity 
is v = 3.5, in line with the experimental results of Renier et al.[6]. 
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4 SURFACE WAVES 



Dowaikh and Ogden[22] established the secular equation for a surface acoustic 
wave propagating in a principal direction of a deformed incompressible solid 
(see also Hayes and Rivlin[25]). Here, for a wave traveling in the direction of 
elongation xi, with attenuation in the ^-direction, it reads in general as 

?7 3 +?7 2 + (2/3i2+272i-7i2)?7/72i-l = 0, where rj = y/ (r/12 ~ pv 2 )/^- 

(16) 

With 7 12 , P12, and 721 given by Eq.fjHJ), the secular equation Eq. ffTElb is 

f(r), e 2 ) = 7] 3 + rj 2 + (3 + 12^e 2 )r/ -1 = 0. (17) 

Clearly now, rj is of the form 77 = 770 + r]2e 2 plus higher-order terms, and no 
term in e 1 . Here, 770 = 0.2956 is the unique positive real root of the cubic 
rj 3 + ?7 2 + 3?7 — 1 = 0. To find 772, we expand Eq.( |T71) . and conclude that 

~dP*\^ ' 12 

» = ~^hr = " wto) "- o - 92 ° 5 ^ « i8) 

dr] 

It follows from this and Eq. (lTE]) 9 that pv 2 = 712 — ^lVoiVo + 2?72e 2 ), or 

pv 2 = 0.9126/i + (3/i + 0.9126A/4)e + (5.642// + 2.071A + 3.554D)e 2 . (19) 

Here, the numbers 3 and 4 are exact, whilst 0.9126 is xr = 1 — i/q, the 
celebrated root of Rayleigh's cubic [21]: x 3 — 8x 2 + 24x — 16 = 0, for surface 
waves in linear, isotropic, incompressible solids. 

The principal Rayleigh wave propagating in the Xi-direction, with atten- 
uation in the x 3 -direction has the same secular equation, due to the equalities 
7i3 = 712, 723 = 7i3, and /?i 3 = (3 U . 

To find the speed of a principal surface wave propagating in the in- 
direction, with attenuation in the xi-direction (or equivalently, propagating 
in the x 3 -direction, with attenuation in the rci-direction), we swap the roles 
played by 712 and 721 in Eq. (fTE]) . In that case we conclude that 

pv 2 = 0.9126/i-(0.2621/i-0.2281A)e + (2.379/i + 1.255A + 3.554 J D)e 2 . (20) 

Finally, the speed of surface waves propagating on the faces parallel to 
the (x 2 , x 3 )-plane cannot be measured, because those faces are clamped by 
the uni- axial tension device. 
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5 CONCLUDING REMARKS 



There are two different paths to the study of non-linear elastic effects. The 
theory of exact non-linear elasticity aims at describing large deformations, 
with no cap on the magnitude of the mechanical fields. For instance, it 
tries to find a strain energy density W which can describe the mechanics 
governing the full extension of a rubber string, i.e. up to stretches of 200- 
800%. The constitutive parameters of W are determined from experimental 
data through a global minimization procedure, by non-linear curve fitting 
over the whole stretch range. However it can be argued that the search 
for a good model encompassing the whole non-linear mechanical behavior 
of elastomers and rubber-like materials is a scientific chimera: so far, many 
forms of strain energies have been proposed, but none has given an accurate 
description of real- world materials for all deformation fields. By contrast, 
the theory of weakly non-linear elasticity takes a step-by-step approach to 
non-linearity, and determines constitutive parameters one by one, as they 
appear successively at the onset of non-linear effects. This goal is quite 
simple to achieve experimentally. However, the mathematics involved in the 
resolution of increasingly non-linear systems of equations are quite intricate 
and sometimes, unnecessary, because exact results already exist in exact non- 
linear elasticity, from which it is easy to deduce weakly non-linear expansions 
at the required order. 

This paper aimed at uncovering one of these overlaps of the two theories, 
specifically the acousto-elastic effect in incompressible elasticity. For example 
(§111), it is possible to measure the speed v of a bulk wave propagating in an 
elongated block. The squared wave speed is linearly related to the elongation 
e, as long as e is small. The coefficients of linear correlation give the second- 
and third-order elasticity constants. When e is increased, the v 2 - e relation 
becomes quadratic, and the coefficient of e 2 gives the fourth-order elastic 
constant. Similar results are found using surface waves (§IV). 

These connections are fruitful for isotropic non-linear elastic solids, and 
others can easily be established in, for example, the cases of reflected and 
refracted waves, torsional waves, plate waves, etc. They provide useful ana- 
lytical benchmarks for the experimental determination of third- and fourth- 
order elastic constants of soft gels and silicones. For biological soft tissues 
however, anisotropy must be included into the analysis, because of the pres- 
ence of collagen fiber bundles. This inclusion can have dramatic consequences 
with respect to small-on- large theory fZE\. It also raises rapidly the number 



10 



of elastic constants to be determined. For instance, the introduction of just 
one family of parallel fibers in an incompressible isotropic matrix raises the 
number of constants from 2 to 7 in third-order elasticity, and from 3 to 13 
in fourth-order elasticity [27]. 

Finally, we mention experimental issues. Obviously, generating, control- 
ling, and measuring a nonlinear (finite-amplitude) shear wave is a delicate 
affair, subject to potential problems arising from beam diffraction, temper- 
ature changes due to wave intensity, microstructural damage, etc. By con- 
trast, infinitesimal (small-amplitude) acoustic waves are routinely produced 
by bulk (BAW) or surface (SAW) acoustic wave transducers. This paper 
shows that the coefficient of nonlinearity for shear waves can be obtained 
thanks to the coupling of an infinitesimal BAW or a SAW with the pre-strain. 
Acousto-elasticity itself is a relatively old branch of physical acoustics, that 
can be dated back to the experiments of Hughes and Kelly[28] in 1953, and 
is now integrated in standard Handbooks, see for example Refs.[29l [30]. It 
is employed to determine experimentally third-order elasticity constants of 
ordinary elastic materials within their elastic limit. Fourth-order constants 
cannot be measured in general for metals, crystals, rocks, and other "stiff" 
elastic solids, because of plastic yield, thermal effects, phase transformation, 
etc. For "soft" elastic solids, however, the elastic limit is far beyond 10% (up 
to hundreds of percent for rubbers). It follows that if the acousto-elastic effect 
can be measured in these solids for infinitesimal pre-strains, then it can also 
be detected, in exactly the same experimental conditions, for slightly larger 
pre-strains. As Figj2] suggests, a pre-strain of a couple percent is sufficient 
to pick up the fourth-order elasticity constant. 
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